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Abstract 


This  report  summarizes  the  research  project  “Theoretical  Investigation  of  3-D  Shock 
Wave  -  Turbulent  Boundary  Layer  Interactions”  during  the  period  1  October  1992 
to  31  July  1996.  The  principal  objectives  of  the  research  are  the  investigation  of  3-D 
shock  wave-turbulent  boundary  layer  interactions  (“3-D  turbulent  interactions”)  for 
complex  configurations  ( e.g .,  asymmetric  shock),  and  the  investigation  of  the  accuracy 
of  Chien’s  k  —  e  turbulence  model,  a  newly  developed  Reynolds  Stress  Equation 
turbulence  model  and  a  k  —  e  model  with  a  new  low  Reynolds  number  correction 
for  3-D  turbulent  interactions,  with  specific  focus  on  the  ability  to  predict  surface 
heat  transfer  and  skin  friction.  Accomplishments  during  the  project  include  the 
development  of  a  Reynolds  Stress  Equation  (RSE)  model  (including  determination  of 
all  constants),  and  computation  of  3-D  asymmetric  and  symmetric  shock  interactions 
at  Mach  4  using  Chien’s  k  —  e  turbulence  model,  the  RSE  model  and  the  k  —  e  model 
with  the  new  low  Reynolds  number  correction.  The  report  also  provides  a  list  of 
publications,  scientific  interactions  and  personnel. 
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Chapter  1 
Overview 


1.1  Introduction 

The  design  of  high  speed  aircraft  depends  critically  upon  knowledge  of  turbulence. 
Aerodynamic  forces  and  thermal  loads  are  strongly  affected  by  the  turbulent  flow 
encompassing  high  speed  aircraft  and  their  subsystems  ( e.g engines).  Presently,  high 
speed  aircraft  design  is  limited,  in  part,  by  insufficient  understanding  of  compressible 
turbulent  flows,  and  inadequate  models  for  predicting  their  behavior  [1,  2]. 

The  study  of  3-D  shock  wave-turbulent  boundary  layer  interactions  (re¬ 
ferred  to  as  “3-D  Turbulent  Interactions”)  is  an  important  area  of  current  research 
in  compressible  turbulence.  Commonly  observed  in  high  speed  flight,  3-D  turbu¬ 
lent  interactions  directly  affect  the  performance  of  high  speed  aircraft.  For  example, 
shock  wave-turbulent  boundary  layer  interaction  is  an  important  consideration  in  all 
supersonic  and  hypersonic  inlet  designs  [3].  Good  understanding  and  prediction  of 
the  shock  wave/turbulent  boundary  layer  interaction  will  provide  important  informa¬ 
tion  to  aircraft  designers.  Knight  recently  reviewed  the  status  of  the  CFD  research 
on  three-dimensional  shock  wave/turbulent  boundary  layer  interaction  for  a  family  of 
geometries  [4].  It  is  now  generally  understood  that  the  shock  wave  structure  and  pres¬ 
sure  field  can  be  accurately  predicted  in  general  for  simple  3-D  shock  wave/turbulent 
boundary  layer  interactions,  e.g.,  sharp  fin  and  swept  compression  corner.  However, 
the  heat  transfer  prediction  remains  a  challenge.  Typically,  predicted  heat  transfer 
shows  a  large  deviation  from  the  experiment.  Such  disagreement  is  believed  to  be 
due  to  the  inadequacy  of  the  turbulence  modeling.  To  investigate  the  effect  of  the 
turbulence  modeling  on  the  simulation  accuracy  of  the  3-D  Turbulent  Interactions, 
three  different  turbulence  models  were  used  in  this  project.  Among  the  turbulence 
models,  the  first  is  the  popularly  used  Chien’s  k  —  e  model  with  low  Reynolds  number 
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correction  [5];  the  second  is  a  RSE  turbulence  model  developed  in  this  project  [6];  and 
the  third  is  the  k  —  t  model  with  the  new  low  Reynolds  number  correction  developed 
under  separate  sponsorship  [7]. 

The  computation  of  the  3-D  Turbulent  Interactions  using  the  RSE  model  in  this 
project  [8,  9]  is  aimed  to  overcome  the  disadvantages  of  the  eddy  viscosity  concept 
inherent  in  algebraic  and  k  —  e  models.  Reynolds  Stress  Equation  (RSE)  turbulence 
models  are  promising  since  they  naturally  incorporate  the  non-local  and  history  effects 
of  the  Reynolds  stress  development.  However,  applications  of  RSE  models  are  mainly 
in  the  research  field  of  incompressible  flow.  The  computation  is  this  project  is  the  first 
application  of  RSE  model  to  a  3-D  shock  wave/turbulent  boundary  layer  interaction 
flow. 

The  two  equation  k—e  model  is  the  still  most  common  choice  since  it  can  in  principle 
better  predict  complex  flowfields  than  algebraic  models  and  is  significantly  simpler 
than  sophisticated  higher  order  closures.  In  this  research  program,  a  k  —  e  model 
with  a  new  low  Reynolds  number  correction  [7]  was  employed  .  This  low  Reynolds 
number  model  was  developed  on  the  basis  of  three  principles,  namely,  1)  the  model 
employs  the  physical  dissipation  rate  e,  2)  the  normal  distance  y  is  avoided,  and  3)  the 
minimum  number  of  modifications  are  introduced,  as  described  by  Speziale  [10].  The 
low  Reynolds  number  modifications  are  1)  incorporation  of  molecular  diffusion  of  k 
and  e,  2)  modification  of  the  turbulent  eddy  viscosity  fix  to  provide  proper  asymptotic 
behavior  near  the  wall,  and  3)  modification  of  the  dissipation  of  e  to  avoid  singularities 
in  the  e  equation  near  the  wall. 

Since  the  structure  of  shock  wave-turbulent  boundary  layer  interactions  for  simple 
geometries  is  generally  well  understood,  attention  has  focused  in  this  project  on  more 
complex  3-D  shock  wave-turbulent  boundary  layer  interaction,  including  the  crossing 
shock  (double  fin)  interaction  (Fig.  1.1). 

This  configuration  has  been  proposed  as  a  possible  hypersonic  inlet  [11,  12].  Most 
of  the  research  work  to  date  has  focused  on  the  symmetric  crossing  shock  interaction 
[13, 14, 15,  16, 17, 18,  19,  20].  Collaborative  experimental  and  computational  research 
has  elucidated  the  wave  and  streamline  structure.  Good  agreement  has  been  observed 
between  computed  and  experimental  surface  pressure  and  flowfield  profiles;  however, 
computed  heat  transfer  has  been  found  to  typically  overestimate  the  experimental 
data  [4,  14,  21].  Various  modifications  to  two  equation  turbulence  models  have  been 
tested  to  improve  the  prediction  of  heat  transfer  [20]. 

There  are  only  a  few  studies  of  the  asymmetric  crossing  shock  interaction.  Garrison 
and  Settles  [22]  [23]  obtained  surface  flow  visualization  and  planar  Laser  Sheet  im¬ 
ages  for  a  serious  of  crossing  shock  configurations  at  Mach  3.9.  Zheltovodov  et  al  [23] 
obtained  surface  pressure,  heat  transfer  and  flow  visualization  for  several  asymmetric 
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crossing  shock  interactions  at  Mach  3  and  4.  The  first  computation  of  the  3D  asym¬ 
metric  crossing  shock  wave/ turbulent  boundary  layer  interaction  was  done  in  this 
research  program  by  using  Chien’s  k  —  e  turbulence  model  [24].  Good  agreement  was 
observed  between  computed  and  experimental  surface  pressure  and  flow  visualization; 
however,  surface  heat  transfer  was  over  predicted  by  100%  within  the  3-D  interaction 
region.  The  newly  developed  RSE  and  k  —  e  model  with  low  Reynolds  number  cor¬ 
rection  were  also  used  to  compute  the  asymmetric  crossing  shock  interaction  flow  to 
examine  the  capability  of  these  models  to  predict  such  complex  flows. 


1.2  Objectives 

The  objectives  of  the  research  program  [25]  are: 

1.  To  study  3-D  turbulent  interactions  in  complex  configurations  including  the 
asymmetric  crossing  shock. 

2.  To  investigate  the  accuracy  of  the  Reynolds  Stress  Equation  turbulence  model 
developed  in  this  program  for  simulation  of  3-D  turbulent  interactions  with 
specific  focus  on  prediction  of  surface  heat  transfer. 

3.  To  investigate  the  accuracy  of  the  k  —  e  turbulence  model  with  a  new  low 
Reynolds  number  correction  for  simulation  of  3-D  interactions  with  specific 
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focus  on  prediction  of  surface  heat  transfer  and  skin  friction. 

The  three  year  research  program  was  initiated  in  October  1992.  In  Chapter  2, 
the  research  accomplishments  for  the  period  October  1992  through  July  1996  are 
presented.  In  Chapter  3,  publications  and  personnel  are  listed. 


Chapter  2 

Research  Accomplishments: 
Oct  ’92  -  July  ’96 


2.1  Turbulence  Models 

2.1.1  Governing  Equations 

The  governing  equations  are  the  Reynolds  averaged  Navier-  Stokes  equations  with 
the  additional  equations  for  the  Reynolds  stresses,  the  turbulent  dissipation  rate  and 
heat  flux  terms.  The  Reynolds-averaged  Navier- Stokes  equations  for  conservation  of 
mass,  momentum  and  energy  are, 


dtp  4  dkpuk  =  0 

(2.1) 

• 

dtpui  +  dkpUiUk  =  -dip  +  dk  (- pu-ul  4  fik) 

(2.2) 

• 

dtpe  +  dk  (pe  +  p)uk-  dk  (- CppT'u'l  -  qk ) 

+dk  (- pu”u'lui  +  uffifcj  4-  dk  (^-^pu'ju'jul  +  u-Tik'j 

(2.3) 

where  dt  =  d  j  dt,  dk  =  d/dxk  and  the  Einstein  summation  convention  is  employed. 
The  overbar  denotes  a  conventional  Reynolds  average,  while  the  tilde  is  used  to  denote 
the  Favre  mass  average.  A  double  superscript  represents  fluctuations  with  respect 
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to  the  Favre  average,  while  a  single  superscript  ’  stands  for  fluctuations  with  respect 
to  the  Reynolds  average. 

In  above  equations,  p  is  the  mean  density,  «,•  is  the  mass-averaged  Cartesian  ve¬ 
locity  with  (ui,v,2,U3)  =  ( u,v,w ),  p  is  the  mean  pressure,  and  e  is  the  mass-averaged 
total  energy  per  unit  mass.  The  following  relations  are  employed  to  evaluate  p  and  e: 

p  =  pRT  (2.4) 

e  =  CvT  +  \-UiUi  +  k  (2.5) 

where  k  is  the  mass-averaged  turbulence  kinetic  energy 

pk  =  ^pu'-u-  (2.6) 

The  mean  molecular  viscous  stress  f,j  is 

2 

Tij  =  dkUkSij  +  p  (diUj  +  djPi)  (2.7) 

where  p  =  p(T)  is  determined  by  the  Sutherland  law.  Similarly,  the  molecular  heat 
flux  is 

qk  =  -TT&r  (2.8) 

rr 

where  Pr  is  the  molecular  Prandtl  number. 

In  the  energy  equation,  the  triple  velocity  correlation  | pu'jUju'l  and  the  velocity- 
molecular  shear  correlation  u^Tik  are  considered  small  and  are  neglected  [26]. 

To  close  the  above  equations,  a  turbulence  model  needs  to  be  introduced  to  deter¬ 
mine  the  Reynolds  stress  —pu'-u'j  and  turbulent  heat  flux  —cppT"u”.  Chien’s  k  —  e 
turbulence  model,  a  RSE  model  and  the  k  —  e  turbulence  model  with  the  new  low 
Reynolds  number  correction  developed  by  Knight  were  used  to  compute  the  above 
Reynolds  stress  terms.  The  following  section  presents  each  turbulence  model  used  in 
this  program  in  detail.  The  above  governing  equations  were  solved  by  using  Roe’s 
upwind  scheme  for  convective  and  pressure  terms  and  central  differencing  for  the  vis¬ 
cous  terms.  Implicit  approximation  factorization  method  was  used  to  integrate  the 
governing  equations  in  temporal  direction  [24,  27,  28]. 


2.1.2  A  Reynolds  Stress  Equation  Turbulence  Model 

A  full  Reynolds  Stress  Equation  (RSE)  turbulence  model  developed  by  Knight [6] 
is  employed.  The  model  is  an  extension  of  an  incompressible  flow  RSE  model  to 
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its  compressible  counterparts  with  some  modification  regarding  the  compressibility 
effects.  The  model  is  presented  below. 

The  equation  for  the  Reynolds  stress  is, 

dtpu'-u'j  +  dkpu'-u'jUk  =  Aij  +  Bij  +  Cij  +  Dij  (2.9) 

where  the  A,j ,  Cij  and  are  the  Reynolds  stress  production  term,  diffusion 

term,  pressure-rate  of  strain  correlation  and  dissipation  term.  They  are  expressed  as: 

A^  =  -pu'ju'ldkUi  -  pu'-uldkUj  (2-10) 

B^  =  dk {  -pu'-Uju’l  +  UjTik  +  u'-Tjk  -pu'-fijk  - puj8ik }  (2.11) 


Cij  =  P  (dju”  +  diu'j)  (2.12) 

Dij  =  -TikdkUj  -  Tjkdku”  (2.13) 

The  production  term  of  the  Reynolds  stress  requires  no  further  modeling  since  the 
Reynolds  stress  and  mean  velocity  are  dependent  variables. 

In  the  RSE  model,  the  diffusion  term  is  modeled  as: 

Bij  =  dk{  Cdl  ~  {pu”u'mdm  (pujttfc)  +  pu'jU^dm  (pU-%)  +  pUkUmdm  (pu[uj)  }} 

+  dk  {  v  {dk  ( pu'-u'j )  +  dj  (pttfttjj)  +  di  }}  (2.14) 

where  i '  =  p/p  and  Cdl  is  a  constant.  This  model  is  an  extension  of  the  incompressible 
flow  model  of  Launder,  Reece  and  Rodi  [29] 

The  correlation  of  the  instantaneous  pressure  and  fluctuating  rate-of-strain  is  mod¬ 
eled  as: 

Cij  =  —  CPl 1 1 pu'-u'j  -  ^pkSij^j  +  C^pk {djiii  +  diUj}  (2.15) 

where  CPl  and  C^  are  constants.  This  is  an  extension  of  Rotta’s  model  [30]  for 
incompressible  flow. 

An  isotropic  dissipation  model  with  compressibility  effect  is  used  to  determine  the 
dissipation  term: 

2  _ 

Dij  —  —p€8ij 


(2.16) 
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According  to  Sarkar  et  al.  [31]  and  Zeman  [32] 

pe  =  p(e,  +  ec) 

where  ec  =  CkCgM?  and  the  turbulence  Mach  number  is  Mt  =  \/2k/a.  Even  though 
this  compressible  flow  effect  for  dissipation  is  incorporated  in  the  computer  code,  it  is 
actually  not  used  and  the  Ck  =  0  because,  for  wall  bounded  flows,  this  effect  is  small 
and  the  above  modification  is  also  questionable  [33]. 

The  conventional  equation  [34]  is  employed  for  e», 


Cg  e2  (  /  p 

UtP^e  dkp^-a^k  =  Cei~pU^UkdkU i  ^  ^  |  (  ^e3 P  ^ 

where  Cei ,  Ce2  and  are  constants. 

The  turbulent  heat  flux  is  modeled  using  a  gradient  diffusion 

-cp^=cppCh^diT  (2.18) 

where  Ch  is  a  constant. 

Table  2.1  presents  the  values  of  the  closure  constants  of  the  present  RSE  turbulence 
model  [35]. 


+  /^&e8J  (2.17) 

hypothesis 


Table  2.1:  Turbulence  Model  Constants 


Quantity 

Value 

Cdl 

0.086 

CPI 

4.325 

Cn 

0.179 

Cei 

1.01 

C« 

1.80 

C' 3 

0.10 

Ck 

0.0 

Ch 

0.0857 

There  are  two  ways  to  treat  the  flow  field  in  the  vicinity  of  the  walls:  one  is  to 
use  wall  functions  and  the  other  is  to  use  the  low  Reynolds  number  correction  to 
resolve  the  viscous  sublayer.  The  wall  function  method  is  chosen  in  this  paper  to  save 
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computational  time.  The  asymptotic  solutions  of  the  above  governing  equations  in  the 
fully  turbulent  region  of  a  2-D  boundary  layer  are  obtained  in  [36].  The  asymptotic 
solutions  (wall  functions)  are: 


u  1 
2A2 


where 


jl?  +  Vb2  +  4A2  sin  A-^—  (log  y+  +  B'k)  — arcsin  (-^= 


2  +  4A2 . 


.  yu*Pw 

y  =  — — 

Pw 

A2  = 


=  0, 
(2.19) 


B  =  -Prt 


Qw  Ua 


CpTwTw 

B'  =  5.0, 


II 

0 

J— 1 

pw 

p 

-1-B—  +  A2(— )2  =  0 

Uoo  \UooJ 

(2.20) 

P  —  RpwTw 

(2.21) 

kiso  {Tw  — 

'I'vjall'j  4"  (1  kzso)  {Qw  Qwall )  —  0. 

(2.22) 

e  ^^ys=0i 

(2.23) 

n  n  _  n 

pu  U  ~  OL\Tw  =  0 

(2.24) 

pv  v  —  a2rw  =  0 

(2.25) 

pw  w  —  a3rw  =  0 

(2.26) 

pu  V  +  Tw  =  0 

(2.27) 

II 

II 

0 

(2.28) 

where  Prt  =  0.9,  «oo  and  T^  are  reference  velocity  and  temperature.  k;30  =  1 
indicates  a  specified  wall  temperature  and  k;30  =  0  implies  a  specified  wall  heat  flux. 
The  constants  «i,  ,  0:3  and  a4  are  defined  as: 


«i  —  (CPl  +  2) 


3  (C„  -  1 


1-5  CnCn) 


(2.20) 
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CP  i  +  2 


a 


2 

4 


«3  =  a  2 


C>2 

_  C„  Cp, 
a  (a2-CP2a)a2 


(2.30) 

(2.31) 

(2.32) 


where  a  =  |(aii  +  a2  +  <*3)- 

These  wall  functions  are  extended  in  a  straightforward  manner  to  3-D  flow  under 
the  assumption  that  the  turbulent  shear  stress  is  locally  aligned  with  the  direction 
of  the  mean  flow  velocity  at  the  point  of  application  of  the  wall  functions.  The 
correction  to  the  wall  functions  due  to  the  pressure  gradient  [37]  is  omitted.  The 
relative  contribution  of  the  pressure  gradient  terms  to  the  Reynolds  shear  stress  at  the 
point  of  application  of  the  wall  functions  [37]  is  proportional  to  y  and  thus  decreases 
with  decreasing  y.  As  indicated  in  [27],  the  computed  flowfield  was  found  to  be 
insensitive  to  the  y- value  of  the  point  of  application  of  the  wall  functions,  and  thus  it 
is  concluded  that  the  pressure  gradient  correction  to  the  wall  functions  would  have  a 
small  effect  on  the  computed  flowfield. 


2.1.3  k  —  e  Turbulence  Model  with  Chien’s  Low  Reynolds 
Number  Correction 


In  this  research  program,  Chien’s  low  Reynolds  number  k  —  e  turbulence  model  was 
firstly  used  to  compute  the  3-D  asymmetric  crossing  shock/ turbulent  boundary  layer 
interaction.  This  model  is  represented  by  the  two-equation  Jones-Launder  (k—e)  [38] 
turbulence  model  incorporating  the  low  Reynolds  number  modifications  of  Chien  [5] 
and  extension  to  compressible  flow  [39] .  The  standard  values  of  the  turbulence  model 
constants  were  employed  [37]  (Table  2.2). 

The  turbulence  model  equations  in  Cartesian  notation  are 


dtpk  +  dkpkuk  =  Pk  +  dk 


—  p[e-\-  2 ukn  2) 


(2.33) 


dtpe  +  8kpeuk  =  Cx^Pk  +  8k  +  /x)  3fce]  -  p |  (C2/e  +  2ukn~2e~Cin+)  (2.34) 


where  p  is  the  mean  density,  uk  is  the  mass-averaged  velocity,  k  is  the  mass-averaged 
turbulence  kinetic  energy,  e  is  the  mean  rate  of  dissipation  of  pt  is  the  turbulent 
eddy  viscosity,  p,  and  v  are  the  mean  molecular  dynamic  and  kinematic  viscosity, 
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Table  2.2:  Turbulence  Model  Constants 


Constant 

Value 

Constant 

Value 

c. 

0.09 

c4 

0.5 

C1 

1.35 

<?k 

1.0 

c2 

1.80 

1.3 

c3 

0.0115 

Prt 

0.9 

respectively,  a*  =  yrw/pw  is  the  friction  velocity,  dt  =  d/dt,  dk  =  d/dxk,  and  the 
Einstein  summation  notation  is  assumed. 

The  production  of  turbulence  kinetic  energy  is 

Pk  =  P t  {9kUj  +  djUk)  dk^j 
where  the  turbulent  eddy  viscosity  is 

fit  =  C>  j  (l  -  e-C 3n+) 

The  dimensionless  damping  function  /  is 

/  =  1  -  0.22e“(fc2/6l/e)2 

The  distance  n  is  measured  normal  to  the  nearest  surface,  and  the  dimensionless 
distance  n+  =  nu*/v. 


2.1.4  k  —  e  Turbulence  Model  with  Low  Reynolds  Number 
Correction  of  Knight 

The  high  Reynolds  number  form  of  the  two-equation  k  —  e  model  of  Jones  and  Laun¬ 
der  [38],  extended  formally  to  compressible  flow  and  employing  the  compressibility 
correction  of  Sarkar  et  al  [31]  and  Zeman  [32],  was  utilized  to  close  the  Reynolds- 
averaged  equations  (2.1)  to  (2.3).  The  low  Reynolds  number  correction  developed 
by  Knight  [7]  employs  the  physical  dissipation  rate  e,  avoids  the  normal  distance  y 
and  utilizes  the  minimum  number  of  modifications  to  the  k  —  e  model,  as  described 
by  Speziale  [10].  The  modifications  to  the  standard  k—e  model  include  incorporation 
of  molecular  diffusion  into  both  turbulent  transport  equations,  multiplication  of  the 
turbulent  eddy  viscosity  by  a  factor  to  provide  the  correct  asymptotic  behavior  of 
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the  turbulent  stresses  close  to  a  solid  boundary,  and  the  inclusion  of  the  dimensionless 
function  /2  for  the  dissipation  term  in  the  solenoidal  dissipation  transport  equation. 
The  governing  equations  are 


dtpk  +  dipkui  =  —pu'-u'-djUi  —  pi  +  di  ( 


ok  4-  tl,d:h 


(2.35) 


dtpev  +  dipuilv  =  -CeiYPu"u'-djUi  -  Ce2f2pj-  +  <9,  (—  + p)  d£ev  (2.36) 

k  k  \  <Te  /  J 


Here  the  dimensionless  function  /2  is  taken  to  be  /2  =  1  —  e~Ce3^  where  Rt  is  the 
turbulence  Reynolds  number  Rt  =  pk2  /fie.  The  dimensionless  function  is  available 
in  tabular  form  as  a  function  of  i2t,  where  f ^  =  0(y~l)  as  y  — +  0  and  /M  — »•  1  as 
y  -*  oo. 

The  standard  values  for  the  high  Re  number  turbtdence  model  constants  were 
utilized  [37].  The  constant  Cu  was  taken  to  be  zero,  since  the  turbulence  Mach  number 
in  the  present  computations  is  small  and  consequently  the  dilatational  dissipation  can 
be  neglected  [40].  The  determination  of  constant  Ce3  is  described  in  detail  in  [7].  All 
turbulence  constants  are  listed  in  Table  2.3. 


Table  2.3:  Standard  k  —  e  Model  Constants 
Constant  Value 


C, 

0.09 

C« 

1.44 

Ce, 

1.92 

Prt 

0.9 

<Tk 

1.0 

cr£ 

1.3 

C. 3 

0.17 

The  boundary  conditions  for  the  turbulence  variables  at  a  solid  boundary  are 


k 

6V 


=  0 

_  2pw  ( dVk\ 

pw  \  dy  ) 


(2.37) 


where  y  is  the  normal  distance  to  the  boundary  and  the  subscript  w  implies  evaluation 
at  the  wall.  These  boundary  conditions  are  exact. 
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2.2  Crossing-Shock  Interaction 

The  following  sections  only  give  a  part  of  the  representative  results.  Detailed  results 
can  be  seen  in  the  cited  references. 


2.2.1  Computation  of  7°  x  11°  Double  Fin  Flow  Field  Using 
the  RSE  Turbulence  Model 

The  asymmetric  double  fin  channel  was  calculated  using  the  RSE  turbulence  model 
[8,  27].  The  two  fins  have  angles  of  7°  and  11°,  respectively.  The  Mach  number  of  the 
upstream  flow  is  3.95.  The  thin  boundary  layers  on  the  fin  surfaces  were  not  resolved, 
and  slip  boundary  conditions  were  applied  at  the  fin  surfaces.  Previous  studies  of 
the  single  fin  interactions  [41]  have  demonstrated  that  the  bottom  wall  shock  wave  - 
turbulent  boundary  layer  interaction  is  essentially  unaffected  by  the  boundary  layer 
on  the  fin.  Of  course,  sufficiently  far  downstream  of  the  intersection  of  the  crossing 
shocks,  the  shocks  will  interact  with  the  turbulent  boundary  layers  on  the  sidewall 
fin  surfaces.  However,  in  the  present  study,  nearly  all  of  the  experimental  data  is 
obtained  upstream  of  the  sidewall  shock  wave-turbulent  boundary  layer  interactions, 
and  consequently  the  omission  of  the  fin  boundary  layers  does  not  affect  comparison 
with  the  experiment  in  this  regard.  Also,  Bardina  and  Coakley  [21]  observed  that 
the  treatment  of  the  fin  boundary  layers  as  laminar /turbulent  affected  the  flowfield 
in  the  immediate  vicinity  of  the  fin/plate  junction  for  a  Mach  8.3  crossing  shock 
interaction.  However,  in  the  present  study,  experimental  data  is  available  only  in  the 
central  portion  of  the  flow,  and  thus  our  comparison  with  experiment  is  expected  to 
be  unaffected  by  the  use  of  slip  boundary  conditions  on  the  fins. 

A  boundary  layer  equation  code  with  the  same  RSE  turbulence  model  [6]  was  used 
to  generate  the  upstream  variable  profiles.  Table  2.4  presents  the  flow  parameters  for 
the  computations. 

As  the  baseline  case,  the  first  computation  for  this  asymmetric  double  fin  flow  used 
a  mesh  size  equal  to  101  x  69  x  49  in  the  x-,  y-  and  z-directions.  The  y+  of  the  first 
grid  point  adjacent  to  the  wall  was  equal  to  28.96  at  the  inflow  boundary.  A  series 
of  computations  were  carried  out  to  test  the  dependence  of  the  present  RSE  model 
on  the  grid  refinement,  upstream  profiles  and  wall  temperature.  Those  results  were 
compared  with  this  baseline  case  and  indicated  that  the  mesh  resolution  of  the  base 
line  case  was  sufficient  [8,  27].  Table  2.5  provides  the  mesh  size  information  for  the 
baseline  case  and  grid  refinement  cases. 

The  experimental  data  was  obtained  by  Zheltovodov  and  his  colleagues  at  the  Insti- 
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Table  2.4:  Flow  Conditions 


Parameter  Experiment  Computation 


S*oo(m) 

0oo  (m) 

Ckoo 

Ptoo  (MPa) 


3.95 

3.033  x  105 
1.12  x  lO"3 
1.3  x  10"4 
7.9  x  10"4 
1.492 
260.4 


3.95 

3.033  x  105 
1.12  x  lO"3 
1.246  x  10"4 
5.91  x  10"4 
1.492 
260.4 


tute  of  Theoretical  and  Applied  Mechanics,  Russian  Academy  of  Sciences,  Novosibirsk 
[23].  Fig.  2.1  shows  the  experimental  measurement  locations  on  the  bottom  surface 
of  the  channel.  The  location  1  is  the  throat  middle  line  which  is  along  the  streamwise 
direction.  Location  2,  3,  and  4  are  in  spanwise  direction  with  the  streamwise  position 
equal  to  46  mm,  79  mm  and  112  mm.  Fig.  2.2  is  the  pressure  distributions  compared 
with  the  experiment  [40]  and  the  results  using  Chien’s  k  —  e  model  [40]  along  the 
throat  middle  line.  In  Fig.  2.2,  the  computed  pressures  using  both  the  RSE  and 
k  —  e  model  along  the  throat  middle  line  agree  well  with  the  experiment  except  at 
the  very  end  where  there  is  a  strong  shock  wave  reflection  from  the  sidewalls.  Since 
the  side  wall  is  treated  as  inviscid,  the  computed  reflected  sidewall  shock  is  displaced 
downstream  of  the  experimental  shock.  This  is  the  cause  of  the  deviation  between 
the  computed  and  experimental  pressure  for  x  >  140  mm. 

Fig.  2.3  is  the  computed  heat  transfer  coefficient  compared  with  the  experiment 
and  the  computational  results  using  Chien’s  k  —  e  model.  The  heat  transfer  coefficient 
Ch  is  defined  as 

c  _ qw{x,z) _ 

h  _  pooUooCp  (Tw{x,  z)  -  Taw(x ,  z)) 

where  qw(x,z )  =  —nwdT/dy  is  the  wall  heat  flux.  In  the  experiment,  the  Ch  has  a 
measurement  accuracy  uncertainty  of  ±10%  to  15%.  Two  computations  are  needed 
to  determine  Chi  i.e.>  an  isothermal  case  and  an  adiabatic  case.  For  the  isothermal 
case,  the  wall  temperature  is  set  to  265° K  and  the  local  heat  transfer  qw(x,  z)  is 
computed.  For  the  adiabatic  case,  the  zero  wall  temperature  gradient  dT/dy  =  0  is 
used  as  the  boundary  condition  and  the  local  adiabatic  wall  temperature  Taw(x ,  z)  is 
computed. 

The  Ch  measurement  locations  are  the  same  as  those  for  the  pressure.  Along  the 
throat  middle  line  shown  in  Fig.  2.3,  the  computational  results  using  k  —  e  model 
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Table  2.5:  Grid  Sizes  for  Mesh  Refinement  Tests 


Parameter  Baseline  in  x  in  y  in  z 


Nx 

101 

201 

101 

101 

Ny 

69 

69 

137 

69 

Nz 

49 

49 

49 

97 

yt 

29 

29 

20.9 

29 

Ayt 

6.21 

6.21 

2.64 

6.21 

yi/Soo 

.00326 

.00326 

.00237 

.00326 

ymax/  &oo 

15.51 

15.51 

15.51 

15.51 

^ymin/^oo 

0.007 

0.007 

0.003 

0.007 

^ymax/ 

1.16 

1.16 

0.621 

1.16 

Ax/Soc 

0.5 

0.25 

0.5 

0.5 

Azminj  <!>oo 

0.2 

0.2 

0.2 

0.1 

Azmax  /  Sqq 

0.4762 

0.4762 

0.4762 

0.2881 

and  RSE  model  agree  quite  well  with  the  experiment  before  the  intersection  of  the 
shock  waves.  Both  of  the  computational  results  are  substantially  higher  than  the 
experiment  after  the  intersection  of  the  shock  wave  interacting  with  the  turbulent 
boundary  layer.  The  maximum  deviation  of  the  RSE  model  from  the  experiment  is 
80%  and  the  k  —  e  model  is  100%.  The  RSE  model  therefore  has  a  20%  accuracy 
improvement. 

Fig.  2.4  shows  the  mean  streamline  structure  with  the  computed  turbulent  kinetic 
energy  contours  developing  from  the  upstream  turbulent  boundary  layer.  The  turbu¬ 
lent  kinetic  energy  was  normalized  by  the  free  stream  mean  flow  kinetic  energy  u ^ 
and  the  x  coordinates  were  normalized  by  the  upstream  boundary  layer  thickness. 
The  two  bulges  on  the  first  turbulent  kinetic  energy  contours  represent  the  initial 
formation  of  the  individual  vortices  generated  by  the  oblique  shock  wave/turbulent 
boundary  layer  interaction.  These  vortices  converge  to  form  a  counter-rotating  vortex 
pair.  The  shock  wave/turbulent  boundary  layer  interaction  enhances  the  computed 
turbulence  intensity  by  increasing  the  maximum  value  of  k  about  3  times  from  the 
channel  entrance  to  exit.  The  turbulent  region  is  also  enlarged  from  the  boundary 
layer  at  the  entrance  to  about  3.5  times  higher  at  the  exit.  These  results  are  conjec¬ 
tures  and  are  not  verified  by  the  experiment. 

The  CPU  time  per  time  step  for  this  baseline  case  with  the  mesh  size  101  x  69  x 
49  is  113.8s  on  the  Cray  C-90.  Approximately  2000  time  steps  are  needed  to  get 
the  converged  solution.  The  total  CPU  time  for  a  converged  solution  is  therefore 
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Figure  2.3:  Heat  transfer  coefficient  distribution  along  the  throat  middle  line 
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0.0000  0.0030  0.0060 
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Figure  2.4:  Turbulent  kinetic  energy  contours  with  streamlines  at  different  locations 
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approximately  63  hours.  The  computer  code  performance  is  160  MFlops. 

In  conclusion,  the  computation  of  the  flow  field  of  a  7°  x  11°  double  fin  using  a 
full  RSE  model  obtained  good  agreement  with  the  experiment  for  the  pressure  field. 
The  computed  heat  transfer  coefficient  achieved  20%  improvement  compared  with  the 
results  from  the  k  —  e  model  with  Chien’s  low  Reynolds  number  correction.  However, 
the  maximum  heat  transfer  coefficient  still  has  a  80%  deviation  from  the  experiment. 
The  computation  has  resolved  the  counter-rotating  vortex  pair  generated  by  the  cross 
shock  wave/turbulent  boundary  layer  interaction,  and  predicts  an  enhanced  turbu¬ 
lence  intensity  between  the  vortices.  The  computer  code  with  the  RSE  turbulence 
model  is  quite  robust. 


2.2.2  Computation  of  7°  x  11°  Double  Fin  Flow  Field  Using 
Chien’s  k  —  e  Turbulence  Model 

The  same  case  was  computed  using  the  popular  Chien  k  —  e  Turbulence  Model  [40] . 
The  following  table  gives  the  details  of  the  computation  mesh  resolution. 


Table  2.6:  Details  of  Computations 


Ref 

«i 

OL2 

Tw 

Nz 

Nv 

Nz 

Ax/Soo 

Ay  min  f  8  oo 

Aymax  1 8  oo 

Azmin  /  <^oo 

Azmax  / ^oo 

Case  la 

7° 

11° 

I 

101 

68 

49 

0.50 

7.64  X  10~5 

0.51 

0.20 

0.48 

0.16 

Case  lb 

7° 

11° 

A 

101 

68 

49 

0.50 

7.64  x  10-5 

0.51 

0.20 

0.48 

0.19 

LEGEND 


Nx 

number  of  points  in  x 

A \rm» 

rms  grid  spacing  at  wall  in  wall  units 

Ny 

number  of  points  in  y 

I 

Isothermal  wall 

Nz 

number  of  points  in  z 

A 

Adiabatic  wall 

The  computed  surface  skin  friction  lines  and  experimental  surface  flow  visualization 
are  displayed  in  Figs.  2.5  and  2.6,  respectively.  In  comparing  these  results  here  , 
there  are  two  important  points  which  must  be  noted.  First,  the  computed  surface 
skin  friction  lines  are  typically  sensitive  to  the  turbulence  model  employed.  Second, 
the  experimental  surface  flow  visualization  is  a  qualitative  technique  which  has  not 
been  calibrated  to  the  actual  mean  surface  shear  stress  direction.  Nonetheless,  the 
comparison  of  these  quantities  is  helpful  in  understanding  the  flowfield  streamline 
structure. 

The  incident  separation  lines  (lines  of  coalescence  [42])  emanating  from  the  fin 
leading  edges  (1  and  2)  are  clearly  evident  in  the  computation  and  experiment.  These 
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separation  lines  are  associated  with  the  incident  single  fin  interactions.  The  computed 
and  experimental  separation  line  angles1  are  in  close  agreement  (Table  2.7). 

The  incident  separation  lines  coalesce  near  the  center  of  the  region  to  form  a  narrow 
band  of  skin  friction  lines  (3)  offset  to  the  left  side  of  the  channel.  This  line  represents 
the  surface  image  of  the  boundary  between  the  left  and  right  vortices  generated  by 
the  incident  single  fin  interactions.  A  second  line  of  coalescence  forms  alongside  on 
the  right  and  farther  downstream  (4)  associated  with  the  appearance  of  a  narrow 
layer  of  reversed  flow  underneath  the  left  side  of  the  right  vortex  as  indicated  in  [40]. 
A  line  of  divergence  is  also  apparent  near  the  right  fin  (5)  associated  with  the  incident 
single  fin  interaction.  A  similar  line  of  divergence  (unmarked)  is  near  the  left  fin. 


Table  2.7:  Separation  Line  Angles 

Configuration  Reference  fii  /3r 
7°  x  11°  Experiment  22°  28° 

Computation  23°  31° 

LEGEND 

fii  separation  angle  for  left  fin 
separation  angle  for  right  fin 


Figure  2.5:  Computed  skin  friction  lines  (7°  x  11°)1  Left  incident  separation  line 
2  Right  incident  separation  line  3  Left  downstream  coalescence  line  4  Right 
downstream  coalescence  line  5  Line  of  divergence  (similar  line  near  left  fin  ) 


1The  uncertainty  in  determining  the  angle  of  the  separation  line  (beyond  the  inception  region) 
from  the  skin  friction  plot  and  surface  flow  visualization  is  ±1°. 
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Figure  2.6:  Experimental  surface  flow  (7°  x  11°) 


Figure  2.7:  Wall  pressure  on  TML  for  7°  x  11° 
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The  surface  pressure  on  the  Throat  Middle  Line  (TML)  is  shown  in  Fig.  2.7. 
Overall,  good  agreement  is  observed  for  x  <  135  mm,  although  the  computation 
underestimates  the  extent  of  the  upstream  influence,  as  observed  in  previous  studies 
( e.g .,  [15]).  The  computed  pressure  does  not  accurately  predict  the  pressure  rise 
(beginning  at  x  —  145  mm)  associated  with  the  shock  reflection  from  the  7°  fin,  since 
the  computation  omits  the  boundary  layers  on  the  fin  surfaces. 

The  computed  and  experimental  surface  heat  transfer  coefficient  Ch  is  presented 
in  Fig.  2.16.  Downstream  of  the  intersection  of  the  fin-generated  A— shock  structures 
(i.e.,  for  x  >  90  mm  on  the  TML)  the  computations  overpredict  the  heat  transfer 
by  approximately  a  factor  of  two.  This  discrepancy  is  believed  attributable  to  the 
limitations  of  the  turbulence  model,  and  is  consistent  with  previous  results  for  the 
symmetric  crossing  shock  at  Mach  8.3  using  a  similar  k  —  e  model  [43,  44]. 

In  conclusion,  the  computed  surface  pressure  displays  good  agreement  with  ex¬ 
periment.  The  computed  surface  skin  friction  lines  and  experimental  surface  flow 
visualization  display  close  agreement  in  the  location  of  the  initial  separation  lines, 
and  are  in  qualitative  agreement  within  the  crossing  shock  interaction  region.  The 
computed  heat  transfer  is  in  error  by  approximately  a  factor  of  two  within  the  region 
downstream  of  the  intersection  of  the  A— shocks  generated  by  the  fins. 

2.2.3  Computation  of  Double  Fins  Using  k  —  e  Turbulence 
Model  with  Knight’s  Low  Reynolds  Correction 

The  7°  x  7°  and  7°  x  11°  fins  with  incoming  Mach  4  were  computed  using  the  k  —  e 
turbulence  model  with  Knight’s  Low  Reynolds  Number  Correction.  The  7°  x  11°  case 
is  the  same  as  the  one  in  the  last  two  sections.  The  7°  x  7°  fin  has  the  same  incoming 
flow  conditions  as  the  7°  x  11°  fin.  Table  2.8  gives  the  details  of  the  computation 
mesh  resolution. 

The  inflow  profiles  were  generated  with  the  boundary  layer  code  [7] ,  which  utilizes 
the  same  turbulence  model.  The  appropriate  inflow  profile  is  considered  to  be  the 
one  created  by  the  boundary  layer  code  which  matches  the  experimental  value  of 
displacement  thickness. 

For  the  7°  x  7°  and  7°  x  11°  configurations,  the  thin  boundary  layers  on  side  walls 
can  be  neglected  since  the  reflected  shock  waves  intersect  the  sidewalls  near  the  exit 
or  not  at  all  [40,  45]. 
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Table  2.8:  Details  of  Computations 


Ref 

a  l 

OL  2 

BC 

Tw 

Nx 

K 

Nz 

Case  la 

7° 

7° 

S 

I 

101 

79 

49 

Case  lb 

7° 

7° 

S 

A 

101 

79 

49 

Case  2a 

7° 

11° 

s 

I 

101 

81 

49 

Case  2b 

7° 

11° 

s 

A 

101 

81 

49 

Case  2c 

7° 

11° 

s 

I 

202 

81 

49 

Case  2d 

7° 

11° 

s 

I 

101 

162 

49 

Case  2e 

7° 

11° 

s 

I 

101 

81 

98 

Case  2f 

7° 

11° 

s 

I 

101 

81 

49 

Ref 

Ax/Soo 

Ay  min  /  <^oo 

Ay  max/  fi<x> 

AZmin/  ^oo 

AZmax  /  o 

Ay2  |  rms 

Ay 2  latter 

Case  la 

0.5 

2.2  x  10~4 

0.5 

0.2 

0.5 

0.55 

0.52 

Case  lb 

0.5 

2.2  x  10~4 

0.5 

0.2 

0.5 

0.63 

0.59 

Case  2a 

0.5 

2.2  x  10~4 

0.5 

0.2 

0.5 

0.70 

0.62 

Case  2b 

0.5 

2.2  x  lO'4 

0.5 

0.2 

0.5 

0.80 

0.72 

Case  2c 

0.25 

2.2  x  10~4 

0.5 

0.2 

0.5 

0.72 

0.65 

Case  2d 

0.5 

1.1  x  10~4 

0.25 

0.2 

0.5 

0.35 

0.31 

Case  2e 

0.5 

2.2  x  10“4 

0.5 

0.1 

0.25 

0.70 

0.63 

Case  2f 

0.5 

2.2  x  10"4 

0.5 

0.2 

0.5 

0.68 

0.60 

LEGEND 


Nx 

number  of  points  in  x 

I 

Isothermal  wall 

Ny 

number  of  points  in  y 

A 

Adiabatic  wall 

Nz 

number  of  points  in  z 

N 

“No-slip”  sidewalls 

Ay2  |rma 
Ay 2  |aver 

rms  grid  spacing  at  wall  in  wall  units 
average  grid  spacing  at  wall  in  wall  units 

S 

“Slip”  sidewalls 
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7°  x  7° 

The  computed  surface  skin  friction  lines  and  experimental  surface  flow  visualization 
for  the  7°  x  7°  configuration  are  presented  in  Figs.  2.8  and  2.9,  respectively.  The 
separation  lines  1  and  2  originating  from  the  fin  leading  edges  are  apparent  in  the 
computation  and  experiment.  The  computed  and  experimental  separation  line  angles 
agree  within  7%.  The  computed  skin  friction  lines  do  not  intersect,  but,  after  changing 
direction,  slowly  converge  towards  each  other.  Two  weak  divergence  lines  3  and  4 
can  be  found  near  the  fin  surfaces.  The  computational  flow  pattern  is  completely 
symmetric;  however,  the  experimental  results  display  a  slight  asymmetric  behavior 
downstream  in  the  vicinity  of  the  centerline. 


to  4^  o  oo 
o  o  o  o 


Z  (mm) 


Figure  2.8:  Computed  skin  friction  lines  for  7°  x  7°: 

1  Left  incident  separation  line 

2  Right  incident  separation  line 
3,  4  Lines  of  divergence 

5  Downstream  coalescence  line 


The  computed  and  experimental  surface  pressure,  normalized  by  the  freestream 
static  pressure  p0 Q,  are  displayed  in  Fig.  2.10  along  the  Throat  Middle  Line  (TML). 
The  uncertainty  in  the  surface  pressure  measurements  is  ±0.5%.  The  computed 
surface  pressure  displays  excellent  agreement  with  experiment. 
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Figure  2.10:  Wall  pressure  on  TML  for  7°  x  7° 
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Figure  2.12:  Taw  on  TML  for  7°  x  7' 
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The  computed  and  experimental  surface  heat  transfer  coefficient  Ch  is  presented 
in  Fig.  2.11.  The  experimental  uncertainty  for  Ch  is  ±10%  to  ±15%.  Close  agree¬ 
ment  with  the  experiment  is  observed.  On  the  TML,  the  heat  transfer  coefficient  is 
predicted  typically  within  25%  in  the  3-D  interaction  region. 

The  adiabatic  wall  temperature  Taw  is  presented  in  Fig.  2.12  along  the  TML.  The 
experimental  uncertainty  in  Taw  is  less  than  0.2%.  Close  agreement  is  observed.  The 
maximum  difference  between  the  predicted  and  measured  Taw  is  less  than  2%. 

7°  x  11° 

For  7°  x  11°  fin,  the  results  using  the  k  —  e  turbulence  model  with  Knight’s  low 
Reynolds  number  correction  were  compared  with  those  using  RSE  model  and  the  k  —  e 
turbulence  model  with  Chien’s  low  Reynolds  number  correction. 

Figs.  2.13  and  2.14  present  the  computed  surface  skin  friction  lines  and  experi¬ 
mental  surface  flow  visualization  respectively.  It  has  been  previously  noted  [40,  43] 
that  the  computed  surface  skin  friction  lines  are  sensitive  to  the  turbulence  model 
employed.  Comparison  of  current  results  with  [40],  Fig.  6,  shows  general  agreement 
as  well  as  a  number  of  substantially  different  details.  Both  incident  separation  lines 
(lines  of  coalescence  [42])  emanating  from  the  fin  leading  edges  (1  and  2)  are  clearly 
observed  in  Fig.  2.13  in  agreement  with  experimental  results  and  previous  simulations 
of  Knight  et  al  [40].  These  separation  lines  are  associated  with  the  incident  single 
fin  interactions.  The  computed  and  experimental  separation  line  angles,  measured 
relative  to  the  x-axis,  agree  within  9%.  However,  contrary  to  the  computation  with 
the  k—e  Chien  model  [40] ,  the  incident  separation  lines  do  not  coalesce  near  the  center 
of  the  region,  but  rather  continue  further  downstream  almost  in  parallel  until  they 
converge  at  x  fa  110  mm  to  form  a  narrow  band  of  skin  friction  lines  (3),  which  is 
offset  to  the  left  side  of  the  channel.  This  is  denoted  in  [40]  as  the  left  downstream 
coalescence  line,  and  represents  the  surface  image  of  the  boundary  between  the  left 
and  right  vortices  generated  by  the  incident  single  fin  interactions.  The  crossflow 
velocity  vectors  near  the  surface  change  direction  at  3.  Lines  of  divergence  are  also 
apparent  near  the  right  fin  (4)  and  left  fin  (5)  associated  with  the  incident  single  fin 
interaction.  In  a  major  difference  with  the  k  —  e  Chien  model  results,  a  second  line 
of  coalescence  (the  right  downstream  coalescence  line)  is  not  present  in  this  compu¬ 
tation.  Consequently,  the  model  does  not  predict  a  secondary  separation  underneath 
the  left  side  of  the  right  vortex  (see  [40]).  The  difference  is  due  to  deviation  in  the 
predictions  of  the  pressure  distribution  in  the  spanwise  direction  obtained  with  each 
turbulence  model. 

The  computed  and  experimental  surface  pressure  along  the  Throat  Middle  Line  is 
displayed  in  Fig.  2.15.  The  computed  and  experimental  surface  pressure  on  TML  are 
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Figure  2.13:  Computed  skin  friction  lines  for  7°  x  11°: 

1  Left  incident  separation  line 

2  Right  incident  separation  line 

3  Left  downstream  coalescence  line 
4,5  Lines  of  divergence 
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Figure  2.14:  Experimental  surface  flow  for  7°  x  11° 

in  good  agreement  for  x  <  135  mm,  although  the  computation  underestimates  the 
extent  of  the  upstream  influence,  as  observed  in  previous  studies  ( e.g .,  [15]  [40]).  The 
computed  pressure  does  not  accurately  predict  the  pressure  rise  (beginning  at  *  =  145 
mm)  associated  with  the  shock  reflection  from  the  7°  fin,  since  all  of  the  computations 
omit  the  boundary  layers  on  the  fin  surfaces. 

The  computed  and  experimental  surface  heat  transfer  coefficient  along  the  TML 
Ch  is  presented  in  Fig.  2.16.  All  three  turbulence  models  overpredict  the  heat  transfer 
by  approximately  a  factor  of  two  downstream  of  the  intersection  of  the  shocks  (which 
occurs  at  x  =  93.7  mm),  with  a  modest  improvement  in  the  computations  performed 
with  RSE  and  present  models. 

The  overprediction  in  Ch  is  actually  an  overprediction  in  qwi  since  a  series  of  stud¬ 
ies  [45,  46]  have  demonstrated  that  the  computed  qw  is  proportional  to  the  computed 
Tw  —  Taw-  A  possible  explanation  is  that  the  turbulence  models  overestimate  the 
effects  of  the  shock-boundary  layer  interaction  on  the  turbulence  production,  thereby 
generating  excessive  turbulence  kinetic  energy  and  overestimating  the  turbulent  ther¬ 
mal  conductivity.  Further  experiments  (in  particular,  measurements  of  the  turbulence 
statistics  within  the  flowfield)  are  needed  to  assist  in  the  identification  of  the  specific 
weaknesses  in  these  models  and  develop  improved  models. 
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Figure  2.15:  Wall  pressure  on  TML  for  7°  x  11° 


Figure  2.16:  Ch  on  TML  for  7°  x  11 
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The  computed  and  experimental  Taw  distributions  in  streamwise  direction  are  dis¬ 
played  in  Figs.  2.17.  The  results  of  the  present  computation  exhibit  excellent  agree¬ 
ment  with  the  experiment,  and  represent  an  improvement  over  the  predictions  by  the 
k—e  Chien  and  RSE  models  [24,  27]. 

The  computed  and  experimental  Taw  contours  are  shown  in  Figs.  2.18  and  2.19, 
respectively.  The  computations  used  the  k  —  e  model  with  the  low  Reynolds  number 
model  of  Knight.  The  results  are  in  qualitative  agreement.  The  computations  predict 
a  modest  increase  in  Taw  downsteam  of  the  intersection  of  the  incident  shocks  and  in 
the  vicinity  of  the  convergence  of  the  incident  separation  lines.  The  decrease  in  Taw 
within  the  expansion  fan  originating  at  x  =  112  mm  is  also  predicted  qualitatively. 


X  (mm) 


Figure  2.17:  Taw  on  TML  for  7°  x  11° 

In  conclusion,  the  computed  surface  pressure  displays  very  good  agreement  with 
experiment  for  the  7°  x  7°  and  7°  x  11°  configurations.  For  the  7°  x  11°  case,  sim¬ 
ilar  very  good  agreement  is  obtained  by  the  k  —  e  Chien  and  full  RSE  models.  The 
computed  surface  skin  friction  lines  are  in  close  agreement  with  experiment  for  the 
initial  separation  lines,  and  are  in  qualitative  agreement  within  the  crossing  shock 
interaction  region.  However,  for  the  7°  x  11°  case,  the  present  model  does  not  predict 
the  secondary  separation  line.  This  feature  is  predicted  by  the  k—e  Chien  model. 
The  adiabatic  wall  temperature  is  accurately  predicted  for  all  configurations.  For  the 
7°  x  11°  case,  the  model  displays  a  definite  improvement  over  the  k—e  Chien  and 
full  RSE  models.  The  computed  and  experimental  Taw  contours  display  qualitative 
agreement,  including  the  modest  increase  in  Taw  downstream  of  the  intersection  of 
the  incident  shocks  and  in  the  vicinity  of  the  convergence  of  the  incident  separation 
lines,  and  the  decrease  in  Taw  in  the  expansion  fan. 
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Figure  2.18:  Computed  Taw  for  7°  x  11° 


EXPERIMENT  TAW/TAW_  TAW_  =239K 


1  2  3  4  5  6  7  8 

0.968  0.981  0.994  1.007  1.020  1.032  1.045  1.058 
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2.3  Summary 

The  k  —  e  and  Reynolds  Stress  Equation  turbulence  models  demonstrated  the  capa¬ 
bility  for  accurately  predicting  certain  features  of  the  flowfield  for  the  crossing  shock 
interaction,  namely,  the  surface  pressure  and  the  surface  flow  visualization.  The  k—e 
model  of  Knight  provided  an  improved  prediction  of  the  surface  adiabatic  wall  tem¬ 
perature.  However,  all  three  models  failed  to  accurately  predict  surface  heat  transfer 
for  the  stronger  (7°  x  11°)  interaction. 

Future  research  should  focus  on  the  investigation  of  more  accurate  turbulence 
models,  specifically  Large  Eddy  Simulation,  for  the  prediction  of  heat  transfer  in 
shock  wave  turbulent  boundary  layer  interactions. 
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Zha,  G.  and  Knight,  D.,  “Computation  of  3-D  Asymmetric  Crossing  Shock  Wave  / 
Turbulent  Boundary  Layer  Interaction  Using  a  Full  Reynolds  Stress  Equation  Tur¬ 
bulence  Model”,  AIAA  Paper  No.  96-0040,  January  1996. 

Gnedin,  M.,  Knight,  D.,  Zheltovodov,  A.,  Maksimov,  A.,  Shevchenko,  A.,  and  Vor¬ 
ontsov,  S.,  “3-D  Crossing  Shock  Wave- Turbulent  Boundary  Layer  Interaction”,  AIAA 
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Paper  No.  96-2001,  June  1996. 

Gnedin,  M.,  Zha,  G.,  and  Knight,  D.,  “3-D  Asymmetric  and  Symmetric  Crossing 
Shock  Wave  Turbulent  Boundary  Layer  Interactions”,  International  Conference  on 
Methods  of  Aerophysical  Research,  Institute  of  Theoretical  and  Applied  Mechanics, 
Russian  Academy  of  Sciences,  Novosibirsk,  Russia,  September  1996. 


3.3  Videotapes 

Narayanswami,  N.,  Knight,  D.,  and  Horstman,  C.,  “Hypersonic  Crossing  Shock  In¬ 
teraction”,  NASA  Ames  Research  Center,  Oct  1993,  11  mins. 


3.4  Presentations 

“Numerical  Simulation  of  3-D  Shock  Wave- Turbulent  Boundary  Layer  Interactions”, 
/ 

Ecole  Polytechnique,  Montreal,  Quebec,  March  1993. 

‘Review  of  Status  of  Computation  of  3-D  Shock  Wave- Turbulent  Boundary  Layer 
Interactions”,  Institute  of  Mechanics  (Moscow  University),  TsAGI  (Moscow  Region), 
CNIIMASH  (Moscow  Region),  CIAM  (Moscow),  August  1994. 

“Development  of  a  Compressible  Reynolds  Stress  Equation  Turbulence  Model”,  Uni¬ 
versity  of  Illinois  at  Urbana-Champaign,  October  1994. 

“Fundamental  Studies  of  3-D  Shock  Wave  Turbulent  Boundary  Layer  Interaction”, 
Boeing  Aircraft  Company  (WA),  June  1995. 

“Numerical  Simulation  of  Shock  Wave  /  Turbulent  Boundary  Layer  Interactions  Using 
Two- Equation  and  Full  Reynolds  Stress  Equation  Models  of  Turbulence”,  Brown 
University,  March  1996. 


3.5  Visitors  to  Rutgers 

Dr.  Doug  Dwoyer,  NASA  Langley  Research  Center,  January  1993. 

Dr.  C.  C.  Horstman,  NASA  Ames  Research  Center,  November  1992,  April  1993. 

Dr.  Len  Sakell,  Air  Force  Office  of  Scientific  Research,  April  1993,  January  1995. 

Dr.  Alexander  Zheltovodov,  Russian  Academy  of  Sciences,  May  1994,  June  1995, 
June  1996. 
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Dr.  Valerie  Zvegintsev,  Russian  Academy  of  Sciences,  May  1994. 
Dr.  Nicolai  Zakharov,  CIAM,  May  1994. 

Prof.  Vladimir  Voul,  Tupolev  Design  Bureau,  May  1994. 


3.6  List  of  Personnel 


Graduate  Student:  Marianna  Gnedin 
Principal  Investigator:  Professor  Doyle  Knight 
PhD  Awarded:  Marianna  Gnedin,  May  1996. 
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